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Abstract 

Fluid instabilities like Rayleigh- Taylor, Richtmyer-Meshkov and Kelvin-Helmholtz instability 
can occur in a wide range of physical phenomenon from astrophysical context to Inertial Con- 
finement Fusion(ICF). Using Layzer's potential flow model, we derive the analytical expressions 
of growth rate of bubble and spike for ideal magnetized fluid in R-T and R-M cases. In presence 
of transverse magnetic field the R-M and R-T instability are suppressed or enhanced depending 
on the direction of magnetic pressure and hydrodynamic pressure. Again the interface of two 
fluid may oscillate if both the fluids are conducting. However the magnetic field has no effect 
in linear case. 
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14 1. Introduction 

15 At the two fluid interface if heavier fluid is supported by lighter fluid, Rayleigh- Taylor instability 

16 (RTI) can occur. When a shock(storng/weak)is passed through the interface of two fluid the inter- 

17 face will be unstable and Richtmyer-Meshkov instability (RMI) can occur. The nonlinear structure 
is resemble like a bubble(when the lighter fluid pushes across the unperturbed interface into the heavier 
is fluid) and a spike (if the fluids are altered) arise due to this kind of fluid instability. RTI and RMI 

20 can occur from astrophysical situation to Inertial Confinement Fusion (ICF). 

21 In ICF, high density fuel is compressed and accelerated towards the origin of the target sphere 

22 by multi KJ laser shock. During shock passage, the interface become unstable which inhibits the 

23 compression in the fusion process. Researchers are searching ways to stabilize these fluid instabilities. 

24 Using Layzer's approximations several authors[l-3] derive the velocity of bubble and spike in linear 

25 and nonlinear domain. Magnetic fields can also be generated due to the ponderomotive force when 

26 the fluids are ionized [4,5]. The effect of magnetic field on Rayleigh - Taylor instability has been 

27 studied in details previously by Chandrasekhar [6]. The growth rate has been found to be lowered 

28 both for continuously accelerated (RTI) and impulsively accelerated (RMI) two fluid interface when 

29 k has a component parallel to the magnetic field [7,8]. 

30 Our paper is addressed to the problem of the temporal development of the nonlinear interfacial 

31 structure caused by RM and RT instability in presence of a magnetic field parallel to the surface of 

32 separation and perpendicular to the acceleration of the two fluids. The wave vector is assumed to lie 

33 in the same plane and perpendicular to the magnetic field. In this type of geometrical situation, no 
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34 effect of the magnetic field is found in the linear approximation. However, in the nonlinear regime, 

35 the effect of magnetic field is predominant. It has been seen that the nonlinear growth rate may be 

36 enhanced or depressed according as the magnetic pressure contribution is either positive or negative. 

37 We have studied analytically and numerically the non linear behavior of the fluid interface in presence 

38 of magnetic field. 

39 2. Basic equations and geometry of the problem 

We have considered two infinite fluids of different constant density separated at y — 0. The heavier 
fluid of density ph is along +ve y axis where as the lighter fluid of density pi lies along -ve y axis. 
Gravity also acts along-ve y direction. The magnetic field acts along z direction. So the Maxwell 

— * — * 

equation is easily valid i.e. V.B = 0. 

After perturbation the finger like nonlinear surface is assumed to be parabolic 



40 Again we consider constant fluid density and hence equation of continuity gives V.v = 0, which 

41 satisfies the irrotational fluid motion. Since we are interested the motion the tip of the bubble, we 

42 can neglect the higher order term of x' l (i > 3) [9]. 





where 



r] > and r] 2 < for bubble 



r/o < and r] 2 > spike 
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Now let us consider the potential function for heavier and lighter fluid, respectively, 
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4> h (x, y, t) = ai (t) cos {kx)e~ k{y - mit)) ; y>0 



(2) 



(x, y, t) = b (t)y + bi(t) cos (kx)e 



k(y-vo(t)). 



(3) 



44 The fluid motion is governed by the ideal MHD equations 



P 



- + (v.V)v 



= -Vp-pg + (JxB) 



(4) 



45 where v = — V0 and J 



VxB 



OB -> -> 

and magnetic induction equation = V x [v x B] 



46 According to our magnetic field consideration 



-(V xB)xB = -(B.V)B - —V(B 2 ) 



(5) 



(6) 



47 Using above relations and substituting in v in Eq. (4) followed by use of Eq. (6) leads to 

48 Bernoulli's equation for the MHD fluid 



P 



-p-pgy-^B 2 + f(t) 



(7) 



49 For RM instability the gravitation acceleration g should be replaced bjg = Au5(t), where Au is 

50 the jump velocity at the interface and S(t) is the delta function. 



si 3. Kinematical and Dynamical boundary conditions at two 

52 fluid interface 

The kinematical boundary conditions are 

I + M.| = («»). (8) 

( V h)x-^ - Wxfa = (Vh)y - {Vl)y (9) 

The Bernoulli's equations for both fluids are 

ph[ - d ~W + l { ^ h)2] - P^- d ~W + = -\9iPh - pi)v + (PH - pi) 

+(f^-|^)] + A(0-/i(*) (io) 

— * 

Further with the help of Eqs. (1) and (2) and the incompressibility condition V ' .Vh{i) = 0, Eq. 
(11) simplifies to 

d[B h{l) (x,y,t)} 

g- t + {v h (i).V)B h(l) = (11) 

53 The above Eqs. [8- 11] give the temporal development bubble at the two fluid interface. 



54 4. Equation for the structure and instability parameters 

55 Substituting the ^ and expanding the velocity terms in powers of the transverse coordinate of x 

56 and keeping up to a; 2 , we obtain the following equations 
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^ = -1(66 + 1)6 (13) 

60 = - ,J* 2 rr fcO! (14) 

(36 + 1) 

61 = TV TT ai ( 15 ) 



6 = k W, 6 = £3 = k ax/yjkg (16) 

57 Where £1,62 and £3 are nondimensionalized (with respect to the wave length) displacement, curvature 

58 and velocity of the tip of the bubble. 

59 Now we are interested in magnetic field induction equations. We set the magnetic field induction 

60 equation to satisfy Maxwell's relation as follows 



B h (x,y,t) = f3 h0 (t) +p hl (t)cos(kx)e- k ^-^; y > (17) 
6i for heavier fluid and for lighter fluid 



Bfayrf = Ao(t) + /9,i(*)cos(fcx)e fc (™<*»; (18) 



6 



62 Using Eq. (11) and expanding the terms up to x 2 for both magnetic field induction equations and 

63 we get 

Pho(t) + CM*) + Phi(t)kfio) cos (kx)e- k ^-^ - k 2 a^ hl e- 2k ^-^ = (19) 

64 for x° : Pho(t) + (3hi(t) = constant = Bho , say similarly for lighter fluid 



Pio(t) + Pn(t) = constant = B m 



(20) 



65 For x 1 



6B h 



6B h (t) (6 



^£ 3 ; 5B h (t) = 



B 



ha 



(21) 



ee which gives 



5B h (t) =5B h (t = 0)exp 



* (e 

?3 



2 " 



(22) 



e? Similarly for lighter fluid 



8B l (& + §)(& + !) XD ^_AW 



5, 



/o 



(23) 



es Hence 
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6B t (t) = 5B t (t = 0)exp 



(24) 



so that &B h (i) < (> 0) if &B ft (i = 0) > (< 0) and 6B t (t) > 0(< 0) if 5B t (t = 0) > (< 0). 
The magnetic pressure difference at the two fluid interface will be 



—Bf(x,y,t) - —B?(x,y,t) = ( 



R 2 R2 
n h0 n l0 \ _|_ ^2 



Now the Bernoulli's Eq.(10) becomes 



^sB h (tm+h-^6 Bl (m 2 



2' 



x 2 (25) 



Ph[- 



dt 



+ ^(V0,) 2 ] - + i(V0,) 2 ] = -g(p h - pi)v + k?^6B h (t)fa + i)x 



(26) 



70 Now we get the following nonlinear equation which describe temporal development of the tip of the 

71 bubble and the velocity of the bubble 



§ =>6+l)6 

|:^(t) _ (6-|) r 

*s fc (0 ~(6 + §r 

(6 + |) (6 + | 



(27) 



rfr " L>(6,r)(66-1) 1 J £>(&,r) 



(66 - 1) r A;V; 



/W 2 



£>(6,0 L 5 



[r— ^<JS ft (*)(26 + 1) + — UB,(f)(2& " 1)] 
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w/iere, r = ty^; r = ^; L>(£ 2 , r) = 12(1 - r)£ 2 + 4(1 - r)& + (r + 1); 



iV(6, r) = 36(1 - r)H + 12(4 + r)£ 2 + (7 - r); V^ (0 = jB 2 /p h(l)m) 



(28) 



72 V fe (/)is the Alfven velocity in the heavier (lighter) fluid. 

73 5. Asymptotic growth rate 

74 To find out the asymptotic value of growth rate of bubble we set d&fdr = which gives £ 2 = —1/6 

75 and at r — > oo integrating the last equation of the set of Eq.(27), 

76 For RTI when lighter fluid is conducting: 



[(&) 



asympibubble 



77 for spike 



/ 2Akg 



I — A kV 2 
l-2{——)^[8B l {oo)\ bubble 
A g 



(29) 



[(&) 



asj/mpj spike 



I 2Akg 
3(T^4)\ 



1 4- A kV 2 
l-2(——)^[SB l (oc)] spike 
A g 



(30) 



78 and for RMI when the lighter fluid is conducting, the asymptotic growth rate is calculated omitting 

79 the second part of the last Eq. of the set of Eq.(27)and integrating,we get 



[(&) 



asympibubble 



4V 2 {1-A) 
\3(Au) 2 (l + A) 



5Bi(oo) bubb i e cot 



4Vf(l 



A) 



3 (1 + A) 
1 + A) \\3(Au) 2 (l + A) 



5Bi(oo) bubb i e 



r 



(31) 



so for spike 
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[(e 



3) asymp] spike 



4lf (l + A) 
\ 3(A-u) 2 (l - A) 



SBi(oo) spike cot 



3 (1 -A) Wfjl + A) 
I- A) ^ 3(A«) 2 (1 - A) 



5Bi(oo) sp ike \ r 



(32) 



81 6. Results and discussions 

82 We have solved the above set of equations using Runge-Kutta-Fehlberg method to describe the tip 

83 of the bubble and the velocity of the tip of the bubble for different cases. 

84 Case 1 

sb Assuming lighter fluid is conducting Bio ^ and the heavier one nonconducting, the hydro- 

se dynamic pressure driven force is suppressed by magnetic pressure. In case of weak shock the RMI 

s? also suppress. For density ratio Ph/pi = r = 1.5, it has been seen that RT instability is suppressed 

ss (Fig. 1). The RM instability is also suppressed in such situation (Fig.2).If the density ratio is increased, 

89 the Atwood number increases, consequently the growth rate increases. The growth rate may decrease 

90 if the density ratio is decreased. 

91 Case 2 

92 If the heavier fluid is conducting and lighter one nonconducting, the magnetic pressure acts along 

93 +ve y direction which increases the bubble growth for both case RT(Fig.3) and RM(Fig.4). 

94 Case 3 

95 If both the fluids are conducting the magnetic pressure difference and hydrodynamical pressure 

96 difference act in different direction and also in opposite phase at the interface. Hence the bubble will 

97 shows oscillatory behavior for both cases. For weak shock the oscillation frequency will be increased 
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98 with the Alfven velocity (Fig. 5). 

99 Applications: 

100 Super Nova explosion starts in a white dwarf as a laminar deflagration and RT instability begins 

101 to act. In white dwarf, magnetic field ~ 10 8 G at the surface and RT instability arising during type la 

102 supernova explosion is associated with the strong magnetic field [10]. In the solar corona, magnetic 

103 field exist in a range of few Gauss to kilo Gauss. The lower limit of magnetic field is ~ 10-20 Gauss, 

104 having temperature 2 x 10 6 k[ll]. Our model suggest that RT instability may show oscillatory 

105 stabilization if the magnetic field is greater than 34 gauss in solar corona. 
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Figure 1: Variation of £1, ^2,bubble growth rate £1) and 5Bi with r for V h 2 = 
for RTI . Initial values are £1 = 0.1, £2 = — 0.1, £3 = 0.1 and r = 1.5 
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Figure 2: Variation of £1, ^2,bubble growth rate £1) and 5Bi with r for V h 2 = 
for RMI . Initial values are as in figl. 
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Figure 3: Variation of £1, £2, bubble growth rate £1) and 5Bi with r for VJ 2 = Ofor 
RTI .Initial values are as in figl. 
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Figure 4: Variation of £1, £2, bubble growth rate £s(= £1) and 5Bh with r for VJ 2 = 
for RMI. Initial values are as in figl. 
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Figure 5: Variation of £1, £2, bubble growth rate &(= £i),SBh and SBiwith r. Initial 
values £1 = 0.1,6 = -0.1, £3 = 0.1,6B h = 5B l = 0.1, B h0 = B l0 = 1.5 and r = 1.5 
blue line for RTI and black line for RMI. 
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